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I study the midgap spectrum of the fermion-vortex system in two spatial dimen- 



(N 



C/3 



X 



sions. The existence of bound states, in addition to the zero modes found by Jackiw 
and Rossi, is established. For a singly quantized vortex, I present complete ana- 
lytical solutions in terms of generalized Laguerre polynomials in the opposite limits 



D ' of vanishing and large vortex core size. There is an infinite number of such bound 



states, with a spectrum that is, when squared, given by, respectively, the Coulomb 
potential and the isotropic harmonic oscillator. Possible experimental signatures of 
this spectrum in condensed-matter realizations of the system are pointed out. 
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O . I. INTRODUCTION 

^ • The important role of topology in determining the properties of physical systems, whether 

oo : 

in cosmology, particle physics, or condensed matter physics is by now well appreciated. This 
role is two-fold. On the one hand, topologically non-trivial structures or excitations such 
as instantons and vortices often arise in various physical systems, from quantum fields of 

00 ! 

^ elementary particles to ordered states of matter, and induce fractional quantum numbers 

relative to the smooth background (i.e. vacuum sector) 1|. Therefore, an understanding of 
their properties is required in order to fully understand the systems in question. On the other 
hand, certain physical systems sustain a novel type of order — topological order — that stems 
from their topological properties In two spatial dimensions, these two roles are closely 
related by the fact that topological excitations acquire fractional exchange statistics [3]. The 
canonical example of this interplay is provided by quantum Hall states with their fractional 
excitations Q and a —tat. degeneracy that depends on the topology of the space in 
which electrons move [5|, I6|j. Due to its global character, a non-trivial topology also offers 
a special sort of protection against local perturbations, which could be useful for coherent 
manipulation of quantum states needed for quantum computation [?]]. 
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A system where such phenomena arise is an ordered state of matter in which fermions 
interact with a vortex in the order parameter. If the topology of the vortex allows an odd 
number of midgap states, the fermion charge becomes fractionalized. When the spectrum is 
particle-hole symmetric, an odd number of midgap states will have zero energy, the so-called 
zero modes, and the fractional charge is found to be |. The number of such zero modes is a 
topological invariant determined by index theorem. An interesting question is whether other 
midgap, non-zero-energy states exist. Apart from its intellectual appeal, a positive answer 
can be useful in detecting such structures in systems where they can potentially arise. 

In superconductors and Fermi superffuids the vortex is carried in the pairing gap function. 
The core states of a vortex in an s-wave superconductor were first calculated by Caroli, de 
Gennes, and Matricon [8], later extended to other limits and studied in great detail 
by others. In chiral superfluids and superconductors, the midgap spectrum was derived by 
Kopnin and Saloma 10J and by Volovik ll|, and the number of zero modes was related to 
the vorticity 121 ] . However, these zero modes are not exact: In a p-wave superconductor, 
there is only one exact zero mode for odd vorticity (l3j |. Moreover, since charge is not a 
good quantum number in a superconductor, the arguments for fractionalization of charge 
do not apply here. Nevertheless, the midgap spectrum has important consequences for 
thermodynamics and transport properties of vortices. 

I study the midgap spectrum of a specific fermion-vortex system, where fermion number 
is a good quantum number and time-reversal symmetry is preserved, unlike quantum Hall 
systems where it is manifestly broken by the external magnetic field. It was originally studied 
by Jackiw and Rossi who showed it supports a number, equal to vorticity, of exact zero 
modes, unlike superconducting vortices. Weinberg [l^ showed this number is protected 
by an index theorem. It was also studied [if]] in a lattice formulation, where similar zero 



modes were obtained. Recently, the same system has been shown [17|, [18| to emerge in the 
low-energy sector of (spinless) fermions hopping on either a honeycomb or a 7r-flux square 
lattice, at half-filling, on the background of a vortex in the dimerized pattern of hopping 
amplitudes. Numerical study of the full lattice Hamiltonian shows 18| the number of exact 
zero modes and the fractional charge bound to the vortex survive beyond the low-energy 



theory. These vortices have 



have also been advanced 



3een shown to obey fractional statistics [19|. Similar proposals 
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22] to realize the zero modes in vortices of other possible 



types of order in graphene-based systems. I note that in Refs. 20j, [22] the zero modes are 
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doubled by the valley index of a graphene layer. Since the vortex inevitably mixes the valleys 
(which are separated by a finite lattice momentum in the Brillouin zone) there is an intrinsic 
splitting of these doubled zero modes and they are not exact |22|. The doubling also makes 
the number of zero modes always even so the fermion number is not fractional. Zero modes 



of Ref. 



21 



occur in vortices of the superconducting order parameter. They are doubled by 
the spin degree of freedom and split by the Zeeman coupling. 

In this paper, I derive the full midgap spectrum of the fermion-vortex system for a 
singly quantized vortex in two physically important limits — those of vanishing and large 
vortex core. Since the splitting of doubled zero modes could be treated as a perturbation, 
the midgap spectrum derived in this paper is relevant in such cases as well. The paper is 
organized as follows. In Sec. [Ill I describe the Hamiltonian. In Sec. IIIIl I derive the equations, 
analyze the general asymptotic solutions, and present the complete midgap spectrum in the 
aforementioned limits. I conclude in Sec. IIVI by examining the relevance of the two limits 
and with brief remarks on the relation to experiment. Some details of the solutions are given 
in an Appendix. 



II. FERMION-VORTEX SYSTEM 



The Hamiltonian is given by 

H = p- a + |m|/3e ix75 , (1) 

where p = — iV is the momentum operator and the matrices a = (7071,7072), P = 7o 
furnish a four dimensional Clifford algebra. For concreteness I use the Weyl representation 
in which 70 = a\ <g> 1, 7fe = — icr 2 <g> cr k (k = 1, 2, 3) and 75 = i7o7i7273 = cr 3 <g> 1 with a k 
the Pauli matrices. The scalar mass m = |m|e lx contains a v- vortex located at the origin, 
i.e., the phase x winds by 2ttv around a loop containing the origin, the amplitude \m\ takes 
a constant value, moo, at infinity and vanishes as rriQr v at the origin. We take a radially 
symmetric form, m(r,9) = \m(r)\e , in polar coordinates (r, 9). 

III. MIDGAP SPECTRUM 

We wish to solve the bound-state spectrum of ([T]), 

Hip = Eip, \E\ < moo. (2) 
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The eigenstate is written as ip — U(ipa,4>b) T where 

U = 



P + P_ 



P- 

is unitary and the projections P± = |(1 ± 03). In this basis, the eignevalue equation ([2]) is 
reduced to 



hlp b = Elp a , rflpa = Eljj b 



h 



m it 



7T m 



(3a) 
(3b) 



with Ti = p x + vp y . We may solve these equations by solving one of two squared equations 



E 



or 



rfhipb = E 2 tp bl ifj a = —hip 



E 



We note that txtx^ = ix^ix = p 2 , and 

htf = 

For a r-vortex, 



Iml 2 + tttt^ \ir, m* 



[m, 7P ] \m\ 2 + 7T7T 



7T 



;m1 = _ ie i(i-^ /d|m| + ^H 
\ dr r 



(4) 



(5) 



(6) 



(7) 



For v > 0, I shall solve the set of Eqs. (HI), which includes the Jackiw- Rossi zero modes with 
ipb = 0. Similarly, for v < we need to solve Eqs. (j5J). Since tr(cx 3 /i/it) = we may take 
^a = (/,r//*) T ,r / = ±l, 



r 2,1 12 P 2A * • i(l-v)e f d \ m \ . u l m l\ « n 

p + m — _£/ ) / — i?7e H ; — 1 / =0. 

\ / \ dr r J 

The dependence on 9 can be eliminated using a two-phase ansatz, 

/(r, 9) = e-Wfxir) + e^ d f 2 (r), £ 1+ £ 2 =v-l. 



(8) 



(9) 



Finally we remove the complex phases by writing fa = e 1Qi 0j, 
±1 = C- Thus, 



i = 1,2 such that ie 1 ^ 1 ^ 2 ) 



1 d d t\ 



r— + -± + \m\ z - E z 0x + (rj 



d\m\ v\m 



r dr dr ' r 2 ' 1 1 J T1 ' \ dr 

1 d d 



+ 



r dr dr r 2 



--— r— + ^ + |m| 2 - E 2 ] 2 + C?7 I + Ui = 0. 



d|m| v\m\ 



dr r 



(10a) 
(10b) 
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At r — > oo, Eqs. ( jTUl) are solved by 0i(r) = C^e "\/ m °° £2r . For r -> we have to the 
leading order, <pi = Dir ai with 

(£l - ai)D ir ai + 2(r](v + l)m D 2 r a2+v+1 = 0, (11a) 
(£l - a\)D 2 r a2 + 2Cr)(v + l)m D ir ai+v+1 = 0, (lib) 

which can be solved (to leading order in r) if either a\ = \£\\, a 2 — a% + v + 1, or a 2 = 
\£ 2 \, ai = a 2 + v + 1. In both cases ai, a 2 > guarantees the regularity at the origin. The 
condition for having a solution over the entire range of r is the smooth matching of the 
two asymptotic solutions. For this we need both asymptotic solutions at the origin. This 
matching would also determine the energy levels. 

I present explicit solutions for the case i\ = i\ = i 2 . For v — 1, we have i\ + i 2 = 0, 
so this condition is satisfied by all the solutions. For v > 1, on the other hand, this means 
h = ^2 = and we find just a small subset of solutions. 

Then, Eqs. (11 01) collapse into one, <pi = —C,r](f) 2 = <p, and take the form of the radial part 
of Schrodinger equation in two dimensions 

1 d d 1,2 



with an effective potential 



~^ r T- + - + V ^ - £ 2 = 0, (12) 
r dr dr r z J 

,9 v\m\ dlml , . 

V(r = m 2 -^^-^^. 13 
r dr 

The sign —^77 is chosen to allow for bound state solutions to exist. 

When v = 1 we may relabel the solutions with I as the angular momentum. From Eq. ([9]) 
we find fi = (e~ l£e e iai + i^e~ 1Ql e l£6 ')0. By taking linear combinations, (ip a )i ± i^e ±2iai (4> a )-e, 
and suppressing overall constant factors the final solutions can be labeled by an angular 
momentum I as 

1 



(14) 



—1 

I now examine the two opposite limits where we can completely solve the spectrum for 
v = l. 



A. Quantum limit 

When the core size of the vortex is vanishingly small, we can take |m| to be a constant, 
moo. Then Eq. ([TBI is reduced to the two-dimensional Coulomb potential, V > (r) = m 2 ^ — 
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vrrioo/r. Taking v — 1, I define 0(r) = r e e ^ r (p(r), with /x = v/mj — -E 2 , and redefine the 
variable z = 2/xr to find 

d 2 , /2£ + l \ d / 1 

d^H~ — V^-^+^-^J^ - (15) 

The solutions to Eq. ([15]) are given in terms of generalized Laguerre polynomials, Ls (z), 
where v = 2£ and s = — £ — \ G N U {0}. So, the spectrum is given by 



tfU(r) = /e-^L^l ( ) , £> = - -, n£N, ( 10) 

2 1 \ 71 J V 71 

where n = 2(s + i) + 1 G N. 

Since s > 0, the maximum value £ can assume for fixed n is found when s — 0, i.e. 
£ < The integrability of |0^ £ | 2 at the origin restricts the minimum value of £ > —^y-- 
To see this, first note that at the origin r\(p^ e (r)\ 2 ~ r q with q = 2£ + 1 + 2p, where p is the 
smallest power of L s . Integrability requires q > —1. For v < 0, p = —uQ(u + s), where 
6 is the Heaviside step function. Thus, q + 1 = £+1- 2£Q(2£ + s) > 0. Since £ + 1 < 0, 
we must have 2£ + s > 0, i.e. £ > — Therefore, \£\ < and the degeneracy of E n is 
# n = n. 

These values of £ also enforce the boundary condition / n ,^(r, + 27r) = (— l) n+1 / n ^(r, 6), 
or conversely, the given boundary condition selects odd or even values of n. The complete 
normalized solutions (including ipj,) are given in the Appendix. 

When v > 1, we find a subset of solutions with 2£ = v — 1, n = 2s + v > v by increments 
of 2, and E> = ±moo-\Jl — In particular, we find only one of the t> zero modes, but we 
also find one non-zero mode per each orbital n. 



B. Core limit 



Deep inside the vortex core we may take the scalar mass to be |m(r)| = m§r v everywhere. 
Then the effective potential, Eq. ( [I3]) . reads 

V<(r) = Ar 2v -Br v -\ (17) 

with A = Wq and B = 2vmo. The spectrum of this potential for general v and E ^ is, to 
the best of my knowledge, not known in closed form. For E = closed-form bound state 



solutions are found [23(, when 



2(v + l)VA 2 v + 1 
again in terms of generalized Laguerre polynomials 

mi ( 2m, 



B ' fcGWU{0}, (18) 



r\%-^ +1 L)r } ( f^r^ ) . (19) 



The bound-state condition, in our case, is v — 1 — 2\£\ = 2(v + l)k, which is satisfied only for 
k = 0. S,nee L<«) = 1, thisjust conn™ the explicit so.ution found originally by Jaekiw and 
Rossi [14{. In contrast to the quantum limit, the non-zero modes are not found for v > 1 in 
this limit. 

When v = 1, the potential V(r) = rn^r 2 — 2mo is that of an isotropic harmonic oscillator. 
So, now the solutions for non-zero E can be found, by shifting B — > B + E 2 in Eq. ffl8|) . as 



< £ = r IV^ 2 L (l^ (mor 2 ); £< = ±v /2m (7i-l), n£N, (20) 

2 

where |£| < n — 1 by increments of 2. Thus, the degeneracy of E n is again g n = n. Since 
t G Z, the boundary condition is necessarily / n ^(r, 6* + 27r) = f n /(r, 9), unlike the quantum 
limit where antiperiodic boundary conditions are also accommodated. 

IV. CONCLUSION 

In sum, I presented, in two physically interesting limits, analytical solutions for the 
midgap spectrum of the vortex-fermion system whose zero modes were found by Jaekiw 



and Rossi 14]. This system has recently emerged as the low-energy theory of certain con- 
densed matter systems, in particular in some ordered states of graphene-based structures. 
I will briefly discuss the relevance of these two limits and remark on possible experimen- 
tal signatures of the midgap spectrum in condensed-matter realizations of the Hamiltonian 
considered here. 

The core size, £, of the vortex can be roughly estimated as £ ~ m^lmQ. The quantum 
and core limits are found, respectively, in the region where r ^> £ and r <^ £. Since the 
zero mode is topologically protected, it is found in both limits. The "high-energy" part, 
where E < m^, is dominated by the quantum-limit spectrum, Eq. (fl6|) . while the non-zero 
"low-energy" part of the midgap spectrum, where E moo, is dominated by the core-limit 
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spectrum, Eq. (|20|) . The low-energy spectrum merges with the high-energy one when the 
orbital quantum number, n = n c ~ ^m^. This implies the low-energy part exists when 
n c > 1, i.e. ^fm^ <^ m^. 

For small enough the low-energy spectrum can be treated as continuous with a 
density of states 

E 3 

n v ' 

By contrast, the density of states of the high-energy spectrum is 

P >(E) = J29nS(E-E>) mlE 



which has a sharp divergence at E — m^. This spectrum must be observable in STM 
experiments that probe the local density of states. It must also have thermodynamic sig- 
natures. For instance, at low temperatures ^m < T C the low-energy spectrum 
renders the specific heat Cy ~ T 4 . Although to detect such thermodynamic dependencies 
experimentally a large enough volume of the system must be occupied by the vortex matter. 
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APPENDIX: NORMALIZED SOLUTIONS IN QUANTUM LIMIT 

Using Eq. (gj) and Eq. dHj) I write 

-ie" ie (£ + f + 4 
Using the properties of generalized Laguerre functions, namely 

A LM(z) = -L<3 L \z), (A.2) 



(A.l) 



dz 

and the recurrence relation 



Lf\z) = L^\z)-L^{z), (A.3) 
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from Eq. (|T6|) we find for level n and 



b)nl 



1 n 



ere ™ 



-ic 



-16 



1 + 



/■ ^ r (l ' +1 Vr , >4./-( I ' +1 Vr1 

!}>(z) — { - z > +L n-i ( z > 



(A.4) 



where Af^e ensures proper normalization, and I remind the reader that v = 2£, s = ! ^ 
and z = ^r. 

n 

The normalization factor can be found by first noting that from Eqs. (jlj) 



i) b \ 2 = —1p* a hh)i) a = |v> a | 



(A.5) 



Using the integral 



z» +l e~ z [LM(z)] 2 dz = (Z/ + S)! (2 g + v + 1), 



We can then find 



nn 



n 
2m r 



£+1 



(A.6) 



(A.7) 



[1] R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 (1976). 

[2] X.-G. Wen, Quantum Field Theory of Many-body Systems (Oxford University Press, 2004). 

[3] M. Oshikawa and T. Senthil, Phys. Rev. Lett. 96, 060601 (2006). 

[4] R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

[5] R. Tao and Y.-S. Wu, Phys. Rev. B 30, 1097 (1984). 

[6] X. G. Wen and Q. Niu, Phys. Rev. B 41, 9377 (1990). 

[7] A. Y. Kitaev, Ann. Phys. 303, 2 (2003). 

[8] C. Caroli, P. G. d. Gennes, and J. Matricon, Phys. Lett. 9, 307 (1964). 

[9] C. Caroli and J. Matricon, Phys. Kondens. Mat. 3, 380 (1965). 
[10] N. B. Kopnin and M. M. Salomaa, Phys. Rev. Lett. 44, 9667 (1991). 
[11] G. E. Volovik, JETP Lett. 70, 609 (1999). 
[12] G. E. Volovik, JETP Lett. 57, 244 (1993). 

[13] V. Gurarie and L. Radzihovsky, Phys. Rev. B 75, 212509 (2007). 
[14] R. Jackiw and P. Rossi, Nucl. Phys. B 190, 681 (1981). 



10 

[15] E. J. Weinberg, Phys. Rev. D 24, 2669 (1981). 

[16] L. F. Cugliandolo, E. Fradkin, and F. A. Schaposnik, Phys. Lett. B 224, 407 (1989). 

[17] C.-Y. Hou, C. Chamon, and C. Mudry, Phys. Rev. Lett. 98, 186809 (2007). 

[18] B. Seradjeh, C. Weeks, and M. Franz, Phys. Rev. B 77, 033104 (2008). 

[19] B. Seradjeh and M. Franz, larXiv:0709.4258l [cond-mat.str-el]. 

[20] I. F. Herbut, Phys. Rev. Lett. 99, 206404 (2007). 

[21] P. Ghaemi and F. Wilczek, arXiv:0709.2626 [cond-mat.supr-con]. 

[22] B. Seradjeh, H. Weber, and M. Franz. larXiv:0806.0849l [cond-mat.str-el]. 

[23] M. Mayrand and L. Vinet, J. Math. Phys. 33, 203 (1992). 



